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A quantum channel can be characterized by a com- 
pletely positive (CP) linear super-operator on the Hilbert 
space of the information carrier. In general this evolution 
is not unitary, so that a pure state loses some coherence 
in its transit through the channel. Various measures of 
a channel's ability to preserve the coherence of its input 
state have been introduced. Arguably, the most useful 
of these are the information capacities 0, 0, El • In con- 
trast, we will focus on the minimum von Neumann en- 
tropy § at the channel output. This quantity is related 
to the minimum amount of noise implicit in the channel: 
when the channel input is a pure state, § quantifies the 
minimum uncertainty occurring in the resulting channel 
output. More precisely, the output entropy associated 
with a pure state measures the entanglement that such a 
state establishes with the environment during the com- 
munication. Because the state of the environment is not 
accessible, this entanglement is responsible for the loss of 
quantum coherence and hence for the injection of noise 
into the channel output. Low values of entanglement, 
i.e., of S, then correspond to low- noise communication. 
Furthermore, the study of S yields important information 
about channel capacities. In particular, an upper bound 
on the classical capacity derives from a lower bound on 
the output entropy of multiple channel uses, see, e.g., 
0, 0|- Finally, the additivity of the minimum entropy 
implies the additivity of the classical capacity and of the 
entanglement of formation [|| . 

Our study of minimum output entropy will be re- 
stricted to bosonic channels in which the electromagnetic 
field, used as the information carrier, interacts with a 
thermal-like noise source. For these channels we analyze 
the following conjecture: that minimum output entropy 
is achieved when the channel input is a coherent state. 
In what follows the rationale and some physical justifica- 
tion for this conjecture are presented. We also consider a 
stronger conjecture — that output states resulting from 
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coherent-state inputs majorize the output states from 
other inputs — which, if true, would imply the minimum 
output entropy conjecture. (Note that the minimum out- 
put entropy problem was previously treated in [jj , which 
reported some of the results that we will discuss.) Ad- 
ditional supporting evidence for the minimum-entropy 
conjecture appears in our companion paper |8j, where 
we show that the integer-order Renyi entropies and the 
Wehrl entropy at the output of the bosonic channel are 
minimized when the channel input is a coherent state. 

In SecUwe present CP-maps for the two bosonic chan- 
nels that will be considered in this paper. The minimum 
output entropy conjecture and its stronger (majoriza- 
tion) version are then stated and explained. In Sec.lTTlwc 
analyze the two channel maps in detail, and develop some 
useful properties of their output entropies. In Sec. lIIII wc 
prove the minimum output entropy conjecture for the re- 
stricted scenario in which only Gaussian states may be 
fed into the channel. In Sec. El we present a collection 
of lower bounds on S. These lower bounds are consistent 
with the minimum output entropy conjecture. More- 
over, in the low and high noise regimes they approach 
asymptotically the upper bounds from which the conjec- 
ture arises. In Sec.[V]we obtain necessary conditions on 
any input state that minimizes the output entropy. We 
demonstrate, in particular, that every coherent-state in- 
put produces an output state which achieves a local min- 
imum of the output entropy. Finally, in Sec. I VII we ad- 
dress the stronger version of the conjecture by exhibiting 
some evidence that output states produced by coherent- 
state inputs majorize all other output states. The paper 
is structured so that Sec. IIIII -Sec. IVII may be read inde- 
pendently. The most technical parts of the derivations 
have been relegated to appendices. 



I. CHANNEL MODELS 

We will consider two bosonic channels — the thermal- 
noise channel and the classical-noise channel — both of 
which belong to the class of Gaussian CP-maps, i.e, 
they evolve Gaussian inputs states into Gaussian out- 
put states ||. Although we will limit our attention to 
single-mode channels, our results can be extended to the 
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multimode case (see, e.g., [To|). 

Thermal-noise channel: Here the signal photons in- 
teract with an environment that is in thermal equilib- 
rium. The channel can be considered to be a beam split- 
ter that couples the input state and the thermal reservoir 
to an output port, with the input-state transmissivity of 
this beam splitter being the channel's quantum efficiency 
n £ [0, 1]. The CP-map £^ for the thermal-noise chan- 
nel is easily obtained by tracing away the thermal state 
of the environment mode (annihilation operator b, with 
average photon number N) from the unitary evolution 



U = cxp 



(bo) — b^a) arctan 



1 — t) 



V 



(1) 



This unitary evolution leads to the beam splitter trans- 
formation 



a — ► U'aU = y/fj a + s/1 — r\ b 
b — ► WbU = y/rjb- VT" 7 ^ a, 



(2) 



where a is the annihilation operator for the channel 
mode. The thermal-noise channel's CP-map is thus 



Tr b [Up® 7* Eft] 



(3) 



where p is the channel input state (of the mode a) and 
n = [N/(N + l)] b1 b /(N + 1) is the environment thermal 
state. The case N — (zero-temperature reservoir) rep- 
resents the pure-loss channel, in which each input pho- 
ton has probability rj of reaching the output. At positive 
temperatures the noise source is active, injecting noise 
photons into the channel mode. For rj = 1, the CP-map 
is the identity: the reservoir does not couple to the out- 
put, hence the 77 = 1 channel is noiseless. 

Classical-noise channel: Here a classical Gaussian 
noise is superimposed on the transmitted field, i.e., the 
classical-noise channel is characterized by the CP-map Q 



Afn(p) 



J d 2 pP n (p)D(p)pD^(p) 



where 



-ImI 2 /™ 



(4) 



(5) 



and D(p) = exp(/.ta^ — p*a) is the displacement operator 
for the mode a. The classical-noise channel's CP-map @ 
is unital — it maps the identity into the identity — whereas 
that for the 77 < 1 thermal-noise channel is not. Never- 
theless, the classical-noise channel can be seen to be a 
limiting case of the thermal-noise channel in which the 
field transformation is replaced by 



(0) 



where p is a classical, complex-valued random variable 
distributed according to Eq. (J5J. The map Af n is then 



obtained from £^ in the limit 77 — i> 1 and N — > 00, with 
(1 - rj)N -> n [see Eq. (JT7J), below]. For n = the CP- 
map 10} becomes the identity, while for n > it injects 
noise photons into the channel. 

A more realistic construct for the classical-noise chan- 
nel is the following. We first propagate the input state 
through a thermal-noise channel with transmissivity n 
and N — 0. Then — to compensate for the propaga- 
tion loss — we employ a phase-insensitive amplifier of gain 
K, = I/77 For 7] = l/(n + 1) the concatenation of 

these two maps exactly yields the classical-noise channel 
Af n (see Eq. (IB51 V 

In response to a vacuum-state input, |0), both the 
thermal-noise and the classical-noise channels produce 
thermal output states given by the density operator 



1 



M 



M + l \ M + 1 



(T) 



with M = (1 — r/)N for the thermal-noise channel £^ , 
and M = n for the classical-noise channel N n - The von 
Neumann entropy of p' is easily found to be S(p' Q ) = 
-Tr^lnp^] = fl (M),withBE3 



g(x) = (1 + a;) ln(l + x) — x lux 



(8) 



for x > and g(0) = (see App.[S}. The output state 
that results from the coherent-state input \a) can be ob- 
tained by the following displacement-operator transfor- 
mation of p' , 



p' a =D{a') p' Dt(a'), 



(9) 



a' = ^ffj a for the thermal-noise channel and 
a' = a for the classical-noise channel [see Eqs. Ijl4|l and 
(TT5|l . below]. Clearly, the states p' a and p' Q have iden- 
tical von Neumann entropies as the former is a unitary 
transformation of the latter . 



A. Minimum output entropy 

We are interested in the minimum von Neumann out- 
put entropy S, which is defined to be 



S(M) 



min S(M(p)) , 

pG/t 



(10) 



where M. is the channel's CP-map and TL is the Hilbert 
space of the information carrier. The concavity of the 
von Neumann entropy 0, implies that § in Eq. {H3J 
can be achieved with a pure-state input. Moreover, given 
two CP-maps, ftA\ and we have that 



S(M 2 oMi) ^ §{M 2 ) , 



(11) 



where the composition (M2 M\)(p) = M.2{M\{p)) is 
the map in which M.2 acts on the output of In- 
equality can be shown by noting that every possible 
input to M.2 on the left-hand side is included in the min- 
imization that is implicit on the right-hand side. 



3 



Conjecture (i): The minimum output entropies 
for the thermal-noise and classical-noise channels are 
achieved by coherent-state inputs, so that 



S = 



g((l-v)N) forS" 



l(n) 



for N n 



(12) 



Discussion: Because the two entropies on the right-hand 
side of (|12|l are achieved by coherent-state inputs, they 
immediately provide upper bounds for the minimum out- 
put entropies of the thermal-noise and classical-noise 
channels, respectively. Thus the conjecture states that 
they are also lower bounds. The conjecture is trivially 
satisfied by the zero-temperature (N = 0) thermal-noise 
channel, because the purity of a coherent state is pre- 
served under the action of the loss map The conjec- 
ture is again trivially satisfied by the n = thermal-noise 
channel, because the map £q sends every input state into 
a thermal state with average photon number N. 

A physical justification for our conjecture resides in 
the fact that in each channel the input state is contam- 
inated by noise from a reservoir whose quantum phase 
is completely random. [The reservoirs are the thermal 
state Tf, for £^ and the classical random source asso- 
ciated with the distribution P n {p) for Af n ]. One thus 
expects that the extraction of any coherence from the 
reservoir — which could be used to reduce the output en- 
tropy below the level when no photons are sent through 
the channel — will be impossible. Some preliminary re- 
sults in this sense were obtained in |l5j . where it was 
shown that, for the thermal-noise map, the linearized en- 
tropy is minimized by a vacuum-state input in the limit 
of low coupling n -C 1 and high temperature N 3> 1. In 
the sections to come we will provide further evidence in 
support of this conjecture. We will also investigate the 
following stronger version of Conjecture (i). 

Conjecture (ii): The output states produced by 
coherent-state inputs majorize all other output states. 
Discussion: By definition, a state p majorizes a state a (a 
property which we denote by p >~ a) if all ordered sums 
of the eigenvalues of p equal or exceed the corresponding 
sums for a pj, i.e., 



p>- a 



(13) 



i=0 



i=0 



where Xi and fii are the eigenvalues of p and a, respec- 
tively, arranged in decreasing order (e.g., A ^ Ai > ■ ■ • ). 
If p >- a then S(p) ^ S(cr), so that Conjecture (ii) im- 
plies Conjecture (i). The converse is generally not true, in 
that S(p) ^ S(cr) does not guarantee that p y a. A nec- 
essary and sufficient condition for a state p to majorize 
a state a is that a be obtainable from p by the action 
of a unital map 0. Conjecture (ii) for the thermal- noise 
channel would then be proved if, for each input state p, 
we could find a unital map C p such that £~ {p) — C p (p' ), 
where p' from Eq. (7J is the output state that is due to a 
vacuum-state input. For the classical- noise channel proof 



we would want a unital map satisfying JV n (p) — C p (p' ). 
We postpone further discussion of Conjecture (ii) until 
Sec. I VII where we will present some evidence that sup- 
ports its validity. The next four sections will concentrate 
on Conjecture (i). 



II. CHANNEL PROPERTIES 

In this section we develop some useful properties of the 
thermal-noise and classical-noise bosonic channels, and 
we provide some insights into their output entropies. 

A. Covariance 

From 0, it is easy to show that the classical- noise 
channel's CP-map is covariant under displacement, i.e., 

N n {D{a) p D\a)) = D(a) N n (p) (a) . (14) 

The thermal-noise channel's CP-map enjoys a similar re- 
lation, viz., 

e%{D{a) P D\a)) = D(^j a) £%{p) D^y/rj a),(15) 

where the ^/rj factor comes from the beam-splitter trans- 
formation (J2J). Moreover, the circular symmetries of the 
probability distribution P n of Eq. JSJ and of the environ- 
ment thermal state tj, imply that both noise maps are 
invariant under the action of a phase shift transforma- 
tion e ^ at<1 . From these two properties plus the unitary 
invariance of the von Neumann entropy 0], it follows 
that any two input states that differ by a displacement 
and/or a phase-shift transformation will produce output 
states with the same entropy. In particular, this means 
that all coherent states produce the same output entropy, 
as discussed previously. 



B. Composition rules 



A complete description of the channels £^ and Af n 
is provided by the transformation of the symmetrically 
ordered characteristic function of the input state 



(16) 



from which p is recovered as Jd 2 /ix(^)_D(— /i)/7r. Evalu- 
ated on the p' , state at the channel output, this function 
becomes 



X'(M) = 



X (VW) e -( 1 -")( JV+1 / 2 )l^l for £ n 



N 



x (p) e-"^ 2 



(17) 



for M n 



As a consequence of the Gaussian character of this evo- 
lution, thermal states (for which i s proportional to 
a zero-mean Gaussian distribution) evolve into thermal 
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FIG. 1: Top panel: graphical representations of the thermal- 
noise and classical-noise channel's CP-maps; the signal pho- 
tons propagate from left to right in these diagrams. Bottom 
panel: composition rules for these channel models; from top 
to bottom are representations of Eqs. 1181 . 11911 . and 12UL 
respectively. 



states. The following composition rules, summarized in 
Fig. ^ follow immediately from Eq. (|17fl : 



Mi 2 oAf nt = M ni +n 2 , 



JVi _ cN' 



in 



(18) 
(19) 



where N' = [772(1 - 771)^1 + (1 - m)N 2 }/{l ~ mm)- 
From Eq. I|18(l we see that concatenating two identi- 
cal classical-noise channels yields another classical-noise 
channel that is twice as random. From Eq. (|T§|> we see 
that concatenating two identical thermal-noise maps re- 
sults in another thermal-noise channel with the same 
reservoir, but whose transmissivity has been squared. It 
is also possible, using Eq. (jT7Jl . to express the thermal- 
noise channel as either a pure-loss channel followed by a 
classical-noise channel, or vice versa: 



£^ — Af(i- V )N £° — N(1-ti)N/ti 



(20) 



From these equations, some useful properties of the out- 
put entropy can be derived. In particular: 

• Because a unital channel increases entropy (the 
output always majorizes the input |lj), Eq. i|18|) 
implies that the entropy S(J\f n (p)) at the output of 
the classical channel is an increasing function of n, 
i.e., for any p and A ^ 



S(Af n+A (p)) > S{Af n (p)) 



(21) 



Because J\f n is unital, we can also infer that the en- 
tropy S(£^ (p)) at the output of the thermal- noise 



channel is an increasing function of N . This fol- 
lows because for any p and A ^ we have that 



S(£^(p)) = 



(22) 



where Eq. ljT%l) and the first equality of Eq. lt2"U|) 

have been used. 

Using (|19[) with Ni = AT 2 = N in conjunction with 
relation (|llf> shows that the minimum output en- 
tropy of the thermal-noise channel is a decreasing 
function of 77, 



for 7/^77 



(23) 



Note, however, that the output entropy S(£^(pj) 
is not a decreasing function of 77 for every p. This 
is because the thermal-noise map does necessarily 
increase the entropy of the input. Consider what 
happens when the channel input is a thermal state 
with average photon number Nq satisfying No > N, 
so that g(N ) > g(N) holds. According to Eq. JT7J), 
the output state is a thermal state with average 
photon number r/No + (1 — 77)^ < No- Its entropy 
is therefore g(nNo + (1 — r))N), which is smaller 
than g(No) and is an increasing function of 77. 

A stronger version of (|23[l can be obtained by us- 
ing relation (|T§|| with Ni ^ N 2 - In this case, (fTT|> 
implies 



>S(£f ) 



for 77' 77 , 



(24) 



and N > 



1-7) 



• The transmissivity inequality in l|23[l can be in- 
verted if the thermal photon numbers are appro- 
priately chosen, viz., 



S(£#') 



for 77 ^ 77' , 



(25) 



where now N' < V-n)"+v'-v 

1 — 7]' 

proven in App. [B] and, together with (1221 
and JI1|), is illustrated in Fig.JSJ 



This relation is 



Using the first equality in Eq. H2U|) , along with (TJ| , 
we can establish the following relation between the 
minimum output entropies of the classical-noise 
map and the thermal-noise map, 



S(f^) = S(Af ( i_„j W o 8°) > §(M(l-r,) N ) 



(26) 



Physically, this says that deleting the pure-loss 
beam-splitter map £ ° can only decrease the output 
noise. An important consequence of (|26() is that if 
Conjecture (i) holds for the classical-noise channel, 
then it must also hold for the thermal-noise chan- 
nel. 
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• The reverse counterpart of l|26|) is given by 

S(Mn) > $(£t7' )/n ') > (27) 

for all n' 6 [0, min{l, n}] (see App.lBlfor the deriva- 
tion). The key consequence of (|27l) is that if Conjec- 
ture (i) is true for the thermal-noise channel, then it 
is must also be valid for the classical-noise channel. 
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FIG. 2: Plot of the minimum entropy regions for the thermal- 
noise channel as constrained by 1221 - 12511 . Each (rj,N) point 
corresponds to a different thermal-noise channel described by 
the CP-map £~ . Given a channel with transmissivity 771 and 
average thermal photon- number Ni, the gray (hatched) re- 
gion represents channels whose minimum entropies are greater 
(less) than S(£^ x ). [The line (r/, 0) has minimum entropy zero, 
and belongs to the hatched region]. In the white regions, the 
composition rules from this section do not establish the rela- 
tion between &(£„ ) and S(f^ x ); in Fig. |5] these regions will 
be partially filled by exploiting the lower bounds that will be 
introduced in Sec. II VI Showing that the upper white region 
is gray and the lower white region is hatched would complete 
the proof of Conjecture (i). Regions 1 and 4 follow from 1221 
and 1231 1 . regions 3 and 6 follow from 1241 . while regions 2 
and 5 are consequences of 1251 . The plot assumes r\\ = 0.7 
and ATi = 0.6. 



and Co = (a, a*). The entropy of pc depends only on its 
covariance matrix 0,0] and is equal to g(VdetT — 1/2). 
Both £^ and Af n transform Gaussian input states into 
Gaussian output states. Moreover, by means of Eq. i|17|l . 
evolution under these CP-maps transforms covariance 
matrices according to 

!r + nl for N n 

(29) 
r,T+ (1 -v)(N + 1/2)1 farf*, 

and first moments according to Co — ► Co f° r Afn> an d 
Co — ► y/vCo f° r ■ The output entropy of a Gaussian 
input state is, hence, equal to g(vdctT' — 1/2), which is 
always greater than or equal to the output entropy of the 
vacuum, i.e., the right-hand-side of Eq. (|12|) . as we now 
will show. 

For the classical map N n it is possible to rewrite l|29|) 

as 

detr' = detr + n(n+ ({Aa,Aa f })) , (30) 

which is always greater than (n + 1/2) 2 because 
({Aa, Aa^}) 1 a nd, from the strong version of the un- 
certainty relation detr ^1/4. In other words, 
we have that 

S{M n { PG )) = g ^det T + n (n + ({Aa, Aat})) _ 1/2^ 

> 9(n) . (31) 

Likewise, we see that Conjecture (i) is true for the 
thermal-noise channel £^ whose input is limited to be a 
Gaussian state because 

detr' = rj 1 detr + (1 - r})(N + l/2)[(l - rj)(N + 1/2) 
+77({A a ,A a t})] , (32) 

which implies detr' ^ [(1 - rj)N + 1/2] 2 . 



III. GAUSSIAN INPUT STATES 



IV. LOWER BOUNDS 



Here we show that Conjecture (i) is true if we restrict 
the channel input to be a Gaussian state pa- Such an 
input has a symmetric characteristic function that is a 
Gaussian form 0, 



x(m) 



exp 



-Co • C f - 2C • r • C f 



C = {ji*,-n), 



which is fully characterized by its first moment Co 
((a), (a')) and its covariance matrix 



r = 



({Aa,Aat})/ 2 ((Aa) 2 ) 
((Aa) 2 ) ({Aa,Aat}>/2 



(28) 



(here ( • ) = Tr[ ■ po\ is expectation with respect to pc, 
Aa = a — (a), and { • , • } denotes the anticommuta- 
tor). The coherent state \a) is Gaussian with T = 1/2 



In this section we present some lower bounds on the 
output entropy. These bounds are consistent with Con- 
jecture (i), and collectively they are asymptotically tight 
in the limits of low and high noise. We will treat the two 
channel models in succession, starting with the classical 
noise case. 



A. Classical-noise channel 

Because g(n) is the output entropy that results when 
the input is a coherent state, g(n) is an upper bound 
on the minimum output entropy S(7V ra ) of the classical- 
noise channel. Four different lower bounds on S(7V" n ) are 
given below. As seen in Fig. |3J bound a is implied by 
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average photon number n — 1, i.e., 



FIG. 3: Bounds on the minimum output entropy of the 
classical-noise channel 7V„ versus average photon number of 
the classical noise, n. Curves a, b, c, and d are the lower 
bounds given in <I351 . J3SJ, id!, and respectively. The 

upper bound u is the function g(n). The minimum entropy 
§(A r n) is constrained by these bounds to lie in the gray re- 
gion, and is required to be an increasing function of n. Con- 
jecture (i) states that S(M n ) = g(n). 



bound d and bound b is implied by bound c. Never- 
theless, we explain all of them because the derivations 
of a and b are simpler. In the limits of low and high 
values of the noise parameter n, it can be shown that 
this collection of bounds is asymptotically tight, i.e., 
lim„^ § c (7V„)/g(n) = lim n _ ) . 00 S d (M n )/g(n) = 1, where 
§>j(Af n ) denotes bound j. 

Lower bound a: By considering the Husimi function 
of the output state, we find that for n ^ 1 



S(M0 ^ g(n - 1) , 



(33) 



by the following argument. Any initial state p can be 
written as [l3 



(34) 



p = d a Q(a) cr(a) , 



where Q(a) = (a\p\a)/7t is that state's Husimi function 
and 

<t(q) = J ^D(X) e >*«-^*+|A| 2 /2 . (35) 
Under the action of the map N n , the state p evolves to 
N n {p) = [d 2 a Q{a) N n (o-{a)) . (36) 



The operator N n (a{a)) is not in general a quantum state. 
However, for n > 1 it is a displaced thermal state with 



A4(a(a)) = D(a) 



D(a)- 



p e 



7t n — 1 

t, 



»</i|Dt(cO 
Dt(a) , (37) 



which has entropy g(n — 1). Lower bound a then can be 
obtained from Eq. I|36|l because Q(a) is a probability dis- 
tribution, and the von Neumann entropy is concave |13| . 

Lower bound b: By considering the Renyi entropy 
S2G0') = — lnTr[(p') 2 ] calculated on the channel output 
p' = Af n {p), we find 

S(.M0 > ln(2n + 1) , (38) 

via two simple steps. As discussed in Q, it is possible to 
show that Sz(p') achieves its minimum ln(2?i + 1) when 
the channel input is a coherent state. Lower bound (|38|l is 
then a trivial consequence of the von Neumann entro py S 
being greater than or equal to the Renyi entropy S2 , |20j . 

Lower bound c: By using a more sophisticated con- 
nection between the von Neumann entropy and the Renyi 
entropy 0, |2^], we find for k ^ 1 an integer and 
n e [(k- l)/2, fc/2] 



S(A/"„) ^ -Afc(n) lnAfe(n) - [1 - A fc (n)] In 
where 



1 - Afc(n) 



X k (n) 



1 - v /l-(A + l)(l-fe/(2n+l)) 



1 



(39) 



(40) 



The derivation of Ij39(l is postponed to App. |0 

Lower bound d: Using the properties of the map Af n 
we find 



S(Af n ) > 1 + bxn, 



(41) 



as we now demonstrate. Consider a pure state \ip) which 
generates an output state with spectral decomposition 



p> =NM){i>\)=Y d lk\ik){ik\, 



(42) 



where {7*:} is a probability distribution and {|7fc)} are 
the orthonormal eigenvectors. From the definition (0J of 
the classical-noise channel's CP-map, we have 



Ik = {lk\p'\lk) 



d 2 pP n (p) ( lk \D(pM . (43) 



The quantity |(7fc|-D(/z)|?/')| is a probability distribution 
over k, and | (7^ \D(p) \ip) | 2 /it is a probability distribution 
over p [23]. Therefore, the convexity of x z for z > 1 
ensures that 



Tr[(p')1 = E ( / d 2 /i J Pn(M)|(7.|^(M)l^)r 

k ^ 

E/ ^ (^M)* |( 7fc |^M| 2 • (44) 
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Because fnO-O i s a Gaussian, it follows that 



,2-1 



(45) 



and the right-hand-side of 14411 can be rewritten in terms 
of the image of \ip) under the action of the map J\f n / Z , 
i.e., 



1 



z w- 



(46) 



This relation can be used to calculate a lower bound for 
the von Neumann entropy by observing that 20] 



S(p') = lim J"W)1 . m ln( Z n^) 

Z->1 Z - 1 Z->1 Z - 1 

= 1 + In n . 



(47) 



Inequality (|47|l applies for any pure state \ip} 7 so we con- 
clude that holds. 





i i i 

E 








W 
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FIG. 4: Bounds on the minimum output entropy of the 
thermal-noise channel £^ as functions of channel transmissiv- 
ity r\ for N = 1/2. Curves A, E and F are the lower bounds 
1481 . 1501 . and 1511 k respectively [Here E is the maximum 
over k of the right-hand side of I5()|l and F is the maximum 
over k of the right-hand side of 15H j. Curves B, C, and Z) 
are the lower bounds lEHJi an d l@D ; respectively, with 

n = (1 — r?)JV. The upper bound u is the function g{{l — rf)N). 
The minimum output entropy S(£^) is constrained to lie in 
the gray region, and is required to be a decreasing function 
of 77. Conjecture (i) states that S(£^) — <?((! — Ti)N). 



B. Thermal- noise channel 

The same techniques that we used to derive lower 
bounds for the classical-noise channel can also be em- 
ployed for the thermal- noise channel £r . The bounds 
we obtain in this case are reported in Figs. 0] [SI and[f)] 



Lower bound A: Repeating the Husimi function cal- 
culation employed above, we find 



S(£%) > g((l - V )N - r,) , 



(48) 



as follows. We replace Eq. Q37[l with 

£?(*(<*)) = (49) 

o ImI 2 

/n // P ( 1 — 77 ) /v — Tf 

which for (1 — 77)^ ^ 77 is a displaced thermal state whose 
entropy equals g((l — r])N — 77). Bound A then follows 
from the concavity of the von Neumann entropy. 

Lower bounds B, C, and D: A simple strategy to de- 
rive bounds on S(£^) is to exploit relation i|26|) . which 
links the minimum output entropy of £^ to that of 
Af(i -r))N - Thus, by replacing n with (1 — 77)^ inequali- 
ties J2BJ), <|39"|) and ijlTJ) of bounds b, c, and d immedi- 
ately become bounds B, C, and D on S(£^), respectively, 
which we have plotted versus channel transmissivity 77 in 
Fig. 0J [Applying this same method to lower bound a, 
we obtain g((l — r/)N — 1), which is not useful, as it is 
already implied by H48[) .] 

Lower bound E: A further lower bound for the 
thermal-noise channel can be derived from the proper- 
ties of the beam splitter; for all integer k it states that 



fc- 
k 



S(C) > 



1 9(^t(1-v)n) 



for 77 < i 



^9{^j[(l-v)N-v+i}) for »H, 

(50) 

the proof appears in App. [DJ Curve E in Fig. 0] is the 
maximum over k of the right-hand side of (|50|l . 

Lower bound F: A more sophisticated version of 
bound E is given by 



> < 



k-l 
k 



' ^ g ((1 - r?)iV) + £S(JV (1 _„ )Ar ) for < 
■g((l-v)N-ri+l) 

+ fe^(^(l— n)N-ri+l/k) 

which is also proven in App. ITT 



for 77 ^ j: 



(51) 

Even though $(Af n ) is 



not known, l|51|l provides a usable lower bound for S(£^) 
when we replace S(J\f n ) with any of the lower bounds 
discussed in Sec. lIV Al Curve Fin Fig.0|is the maximum 
over k of the right-hand side of 1)5 ip. If Conjecture (i) is 
true, then l|51|l becomes an equality for 77 ^ 1/k , i.e., this 
bound is tight. The same is not true for the bound J5UJ. 



V. NECESSARY CONDITIONS FOR MINIMUM 
OUTPUT ENTROPY 

In this section we study the conditions that an input 
state must to satisfy in order to provide a (local or global) 
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FIG. 5: Same as Fig. Q]but for different values of the param- 
eter N: in the left plot N = 0.1, in the right plot N = 10. 
At both high and low average thermal photon numbers, the 
greatest of these lower bounds approaches the upper bound. 




FIG. 6: Plot of the minimum entropy regions for the thermal- 
noise channel for the same parameters used in Fig. |5] Here, 
the region L has been added to the hatched region of Fig. [5] 
by comparing the lower bounds A-F of the channel E^ 1 with 
the upper bound g((l — r/)N) for generic r\ and N: if the 
former is larger than the latter, we can conclude that §(f,^ 1 ) 
is greater than H(£^). Likewise, the region U has been added 
to the gray region of Fig. [5] by comparing the upper bound 
on §(£^ 1 ) with the lower bounds on S(£^). 



minimum for the output entropy. We will show that 
coherent-state inputs satisfy all of these conditions for 
both the classical-noise channel and the thermal-noise 
channel. We begin, in the first subsection, by deriving a 
local minimum condition from the directional derivative 
of the output entropy. Then, in the second subsection, 
we generalize this result into a more stringent condition 
for minimality. 



A. Local minimum condition 

Given a CP-map A4 , a necessary condition for an input 
state <To to provide a local minimum of the output entropy 
S(A4(p)) can be obtained from the directional derivatives 



of this quantity 0. Given t G [0, 1] and a generic state 
cr, this condition requires that 



d_ 

Of 



S(M(a(t))) 



(52) 



Tr 



{[M(<r )-M(a)]\jxM(tT )} 



> 



where a(t) is the mixed state (1 — t)ao+ta. For both the 
thermal-noise map £^ and the classical- noise map Af n , 
this condition is satisfied by the arbitrary coherent state 
(To = l Q: )( Q: |) as can be shown by using covariance prop- 
erties of the noise (see Sec. Ill A|) to rewrite the output 
entropy derivative as follows: 



t=o+ 

Tr[a^a{M{d-)- p' )]bx 
M 



= Tr[(p' -M(a))hxp' ] 
M 



= £Tr[a' a cr] In 



M + 1 



M 
>0, 



1 



(53) 



where ( = I for Af n and £ = rj for £^ , p' Q is the output 
state Eq. (JJJ) generated by using a vacuum state input, 
and a is the state D'(a) a D(a). The last equality in 
(15311 derives from 



Tr [a f a M (a)] = CTr[a f a a] + M 



(54) 



which holds because M is the average number of photons 
in the state p' . Physically, the inequality in (|53|l is a 
consequence of the vacuum's being the input state that 
produces the output state with lowest average photon 
number. 



B. Eigenvalue minimum condition 

In deriving (|52|l we required that the entropy be locally 
increasing when moving along the trajectory a(t), whose 
intermediate states are all mixed. A more stringent re- 
quirement follows from using the pure-state trajectory 
ce = \ag)(ag\ with 



\<ro) 



(55) 



where ctq = |cro)(co| is the putative minimizing state for 
the output entropy, and |oj_) is any state that is orthog- 
onal to | (To)- Expanding the output state entropy in a 
Taylor's series, we have 



S(M(ag)) = 
S(M(a )) 



_0 

00 



7 S(M(a e ,)) 



(56) 



g'=o 



[The entropies of the channels we are considering are dif- 
ferentiable.] If the state |<7o) is a local minimum for the 
output entropy, then the term that is linear in 6 must 
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vanish. This requirement implies the following necessary 
condition for local minimality 



0_ 



S(M(a e )) 



1 = 



-2Rc(a±\M* (lnM(a )) \a ) = 



(57) 



where M* is the dual map associated with 
Ai, such that for any two operators A and B, 
Tt[AM(B)} =Tr[M*(A)B]. Inasmuch as this condition 
must be valid for all \<j±), we can conclude that the 
operator 



T{a ) = -M*(\nM(a )) 



(58) 



has |<7o) as an eigenvector. The properties of T{(Jq) guar- 
antee that the eigenvalue associated with |tro) is its out- 
put entropy. Shor and Ruskai [24j have a different way 
of introducing the operator T{ao) to study the minima 
of a channel's output entropy; their approach does not 
require the output entropy to be diffcrcntiable. 

For the channels we arc considering, the vacuum input 
evolves into the thermal state p' Q of Eq. J7|), so that 



.F(|0)(0|) = M*{(Ja In 



M + 1 

M 



+ ln(M + l)). (59) 



In particular, for the classical-noise channel this quantity 
simplifies appreciably because Af n is its own dual. Indeed, 
because Af n (a^a) = a^a + n, we find 



jF(|0)(0|) = g{n)t + o)a In ■ 



(60) 



which shows that the vacuum is an eigenvector of 
JF(|0)(0|), and hence satisfies the local minimum con- 
dition of Eq. l(57|) . The positivity of the operator a^a 
implies that the vacuum is also the eigenvector with min- 
imum eigenvalue. Using the definition of T, Eq. (|6U[1 can 
be used to express the output entropy of a general pure- 
state input p as 



S(Af n (p)) = g(n)+E In- 



1 



S(Af n (p)\\p' ), (61) 



where E =Tr[a^ap] is the average photon number of the 
input p and S(p\\\p2) =Tr[/>i (lnpi — ln/^)] is the rela- 
tive entropy between states p\ and pi . This equation 
allows us to restate Conjecture (i) in the form 



E\sx- 



> S{M n {p) || p' ) 



(62) 



Proving this relation for all p is equivalent to proving 
Conjecture (i). For coherent-state inputs, we can use the 
covariance properties of the noise under displacements, 
to show that 



F(\a){a\) = D(a)F(\0)(0\)rt(a) , 



(63) 



which guarantees that the coherent state |a) is an eigen- 
vector of T(\a)(a\). If we could prove that the coherent 



states are the only ones which satisfy this condition, we 
would have succeeded in proving the conjecture: the co- 
herent states would be the only states that satisfy the 
necessary condition (|57|) for minimality. Unfortunately, 
such is not the case because Fock states \n) are also 
eigenvectors of the corresponding !F(\n}(n\); this follows 
from the states Af n {\n){n\) being diagonal in the Fock 
basis [25| . Fock states other than the vacuum are not, 
however, minima for the output entropy, as discussed in 
Sec. IVT1 Note that condition was first introduced 
in Ref. [3, where it was claimed that Conjecture (i) was 
proven. A more careful analysis of reveals a funda- 
mental missing link in that proof: even though it is shown 
that number states satisfy condition f5%|l and that the 
vacuum is the number state with the lowest output en- 
tropy, it is not proven that they are the only states that 
satisfy this condition. Hence, there is still a possibility 
that another state (with output entropy lower than the 
vacuum) might satisfy the condition. 

The case of the thermal- noise channel can be treated in 
a similar manner, showing that here too coherent states 
are eigenvectors of the corresponding T operators. In 
this case, the dual map of £ ^ can be written as 

(£fy(p) = Tr b [(l <g> Tb )U\p® t b )U] , (64) 
which is unital and satisfies 

{£*)*{a)a) ^^a+il-^N . (65) 
Using these properties and Eq. I|59|l . we find 



•F(|0)(0|) = .g((l-7?)A0l+7?a T a ln 



(l-rj)N 



(66) 



As in the case of Eq. I|60() , the vacuum is the eigenvector 
with minimum eigenvalue of the operator T . Moreover, 
Eq. H63JI applies here, i.e., the coherent state \a) is an 
eigenvector of ^-"(|a)(ci!|), and Eq. H61|) becomes 

S(£»(p)) = g((l- ?? )JV)+^ln (1 ~^ + 1 

-S(^(p) || p' ) . (67) 



VI. MAJORIZATION 

In the previous sections we gave some justifications in 
support of Conjecture (i). Here we focus on the stronger 
version of this conjecture, i.e., Conjecture (ii). We begin 
by presenting evidence that the output states generated 
by coherent states majorize the ones generated by Fock 
states. 



A. Fock-state inputs 

Because the output states generated by the coherent 
state are all unitarily equivalent, we can focus on the 
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ordered eigenvalue sums for the thermal state p' from 
Eq. 0, i.e., 



M 



+ 1 ^ 

i=0 



A/ 



M + l 



= 1 - 



M 



M+l 



9+1 



(68) 



The ordered partial sums i|68[l for all g must be compared 
with their Fock-state-input counterparts. In the case of 
classical noise, these can be numerically evaluated by ob- 
serving that, for a Fock input state \k), the output is 
diagonal in the Fock basis and takes the form [2j| 




KQk){k\) = X>|i)(i|, 



i=0 



where 



min(fc,i) 

E 

i=o 



,fe+i-2j 



( n + l)fe+i+i 



(69) 



(70) 



[Note that in the case k = 0, Eq. 1)69(1 reduces to the 
vacuum evolution Eq. JJJ]. Evaluation of the ordered 
partial sums is tedious but can be performed analytically. 
In particular, for k = 1 and n 1, the ordered partial 
sums {Sg} of the first q + 1 eigenvalues are 



E, = l- 



1 



n(ra +1)/ \n + 1/ 



(71) 



which, for all q, are smaller than the corresponding sums 
in Eq. iJBSJ. The case n < 1 is analogous: for sufficiently 
large values of q the sum is the same as in Eq. I|71[) . while 
for small q it is given by 



E, = l- 



1 



q + 2 
n(n +1)/ \n + 1 



q+2 



(n + 1) : 



r,(72) 



which again is smaller than the sum in Eq. (|68|) . In Fig.[7| 
a numerical comparison between Eq. (|68l) and the ordered 
partial sums of the eigenvalues from Eq. I)69|) is presented 
for two specific cases; p' majorizes Af n (\k) (k\) in both. 

The same analysis can be repeated in the case of the 
thermal channel, observing that 



£^\k)(k\) = J2Pm H(m\ 



m=0 



where {p m } is the binomial distribution 
k 



Pn 



k — rn 



(73) 



(74) 



Using Eq. (|69l) and the decomposition 120|) the thermal 
evolution of the Fock state \k) can be calculated as 



k 

E 

m=0 



£^{\k){k\)=^p m U [1 _ r])N {\m){m\) 



(75) 



FIG. 7: Majorization analysis of the map J\f„. Plots of E 9 , the 
sum of the q + 1 largest eigenvalues, for the thermal state p' 
(gray staircase) and the evolved Fock state J\f n (\k) (k\) (black 
staircase). In both cases the thermal state majorizes the 
evolved Fock state, which implies that the entropy of the for- 
mer is smaller than the entropy of the latter p|. The insets 
show the eigenvalues {Xi} as functions of the photon number 
i. Here k = 1, n = 0.85 for the left plot, and k = 6, n = 0.85 
for the right plot. 



Note that the output is again diagonal in the Fock ba- 
sis. Moreover, if the output of the vacuum majorizes the 
output of the other Fock states for the classical channel 
J\f n , Eq. I|75|) can be used to prove that this must also be 
true for the thermal channel £^ . In fact, if jV n (|0)(0|) 
majorizes AT n (\m) (m\), then there exists |l| a unital map 
C m such that Af n (\m){m\) = C m (N n (\0)(0\)). Now, be- 
»))A r (l )(°l)' E< F implies that 



cause £,™(|0}(0|) =7V (1 



C(|fc)(fc|) 



5> m C m (^(|0)(0|)) 



(76) 



The convex sum of unital maps is a unital map, hence 
£^(\k)(k\) is obtained from a unital transformation of 
£^(|0)(0|), which implies Q that the latter majorizes 
the former. 



B. Arbitrary input states 



Further insight into Conjecture (ii) is provided by gen- 
eralizing our analysis to the case of an arbitrary pure- 
state input. Because it is sufficient to establish Conjec- 
ture (ii) for the classical-noise channel, we shall only con- 
sider that case. As shown in App.[El see also [2{|, when 
an arbitrary pure state is fed into the classical- noise 
map Af n , the resulting output state p' has the following 
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Fock-state representation, 



C. Simulated annealing optimization 



(k + l\p'\k) = ^j 

„ U + k + i y. 



M 



E 



j\k\ (i + n y+ k + l + 1 

x F(- h -k;-(j + k + l);l-n- 2 ), (77) 

for k, I > 0, where {ip n } are the Fock-state coefficients of 
\ip), and F(a, [3; 7; z) is the hypergeometric function. 

To probe the output-state eigenvalue behavior associ- 
ated with an arbitrary pure-state input with up to 10 
photons, we used the following procedure. Eleven com- 
plex numbers, { </>„ : < n < 10 }, whose real and imagi- 
nary parts were randomly distributed on the interval [-1, 
1] were used to generate a pure state, via 
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|V) = 5> n |n)/ 



n=0 



\ 



10 

Ei 

71=0 



(78) 



Using the {ipn} for this state we diagonalized p' — found 
from Eq. I|77|l truncated to the Fock states { |n) : < n < 
40 } — and then calculated the ordered eigenvalue sum, 
S g , for < q < 40. In all 100 trials of this procedure, 
we found that the output state generated by the vacuum 
majorized the output state produced by the arbitrary 
pure state. Figure [S] shows this comparison for five of 
our 100 trials. These particular input states had average 
amplitudes 1.805 - 0.002i, 1.318 - 0.340i, 1.596 + 0.404i, 
1.810 — 0.255i, 1.546 + 0.276i, and average photon num- 
bers 5.059, 3.375, 4.976, 4.748, 4.163, respectively. 




FIG. 8: Majorization analysis of the map N n for n = 0.85. 
Plots of the sum of the q+1 largest eigenvalues for the thermal 
state po = A/",i(|0)(0j) and the evolved pure state A/"n(|V>) (^l) 
for five randomly-generated \tp). In all cases the thermal state 
majorizes the evolved \4>) state, which implies that the en- 
tropy of the former is smaller than the entropy of the lat- 
ter Q. 



As a final test of Conjecture (ii), we used simulated 
annealing |27| — a well known technique for finding global 
extrema — to minimize the classical-noise channel's out- 
put entropy. As in the previous subsection, the input 
state was truncated to lie in span{ \n) : < n < 10}, 
the output states were constrained to lie in span{|n) : 
< n < 40}, and we limited our consideration to the 
classical-noise channel. A variety of initial pure-state in- 
puts were employed, in conjunction with an exponential 
cooling schedule. In all cases, the resulting minimum 
output entropy was extremely close to that achieved by 
a vacuum-state input. Indeed, in every case the asso- 
ciated input state — at the end of the simulated anneal- 
ing iterations — was very nearly a coherent state. Fig- 
ure OH shows the progression of output entropy values 
for the n = 0.85 classical-noise channel when the sim- 
ulated annealing procedure was initiated with the Fock- 
state input |6) and 400 iterations were performed. The 
initial output entropy in this run was 3.754; the final 
output entropy in this run was 1.846. The latter is 
very close to g(0. 85) = 1.841, which is the output en- 
tropy for a coherent-state input. The final input state 
|'0)4OOj after the 400 iterations, had mean amplitude 
—0.116 + 1.861«, average photon number 3.47, and 99.9% 
overlap with the coherent state \a) for a = — 0.12 + 1. 88i, 
viz., |(a|^)4oo| 2 = 0.999. For this Fock-state input, we 
found that the output-state eigenvalues at every iteration 
majorized those for preceding iterations, see, e.g., Fig. [5] 
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FIG. 9: Left: simulated annealing minimization of the 
classical-noise channel's output entropy for n = 0.85. Plot of 
the output entropy S(Afn(p)) versus iteration number when 
the minimization algorithm is initiated with the Fock-state 
input |6). Right: majorization analysis of the map J\f n for 
n = 0.85. Plots of E 9 (the sum of the q + 1 largest eigen- 
values) for the thermal state p' — 7V n (|0)(0|), the evolved 
Fock state |6), and the output states obtained after 100, 200, 
and 400 iterations of the simulated annealing algorithm. The 
curves for the vacuum-state input and the output state after 
400 iterations are indistinguishable on this scale. 
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VII. CONCLUSIONS 

We conjectured that the minimum entropies at the out- 
put of two Gaussian bosonic channels (with thermal or 
classical noise) are achieved by inputs that are coher- 
ent states. Physically, this conjecture is plausible: the 
complete absence of correlation between the input state 
and the channel's environment state would seem to for- 
bid the existence of an input whose injection reduces the 
output entropy to a level lower than that achieved when 
no photons are transmitted. In support of our conjecture, 
we presented four separate arguments. First, we proved 
that the conjecture is true when we restrict the analy- 
sis to Gaussian-state inputs. Second, we established a 
suite of lower bounds on the minimum output entropy, 
which are all compatible with the conjecture and which 
show that the conjecture is asymptotically correct at low 
and high noise levels. Third, we studied local minimum 
conditions on the output entropy; input coherent states 
were shown to be local minima that satisfy the operator 
identities which are necessary conditions for minimality. 
Fourth, we analyzed a stronger version of the conjecture, 
namely that the output state produced by a coherent- 
state input majorizes all other output states. In sup- 
port of this stronger conjecture, we presented evidence 
for number-state inputs and randomly-selected inputs. 
In a companion paper we show that coherent states 
minimize the output Renyi and Wehrl entropies for the 
classical-noise and thermal-noise channels, thus lending 
further credence to the conjecture in the present work. 



APPENDIX A: VACUUM OUTPUT ENTROPY 

Here we derive the output entropies of the channels J\f n 
and £^ for vacuum-state inputs. Both channels evolve 
the vacuum into the thermal state p' of Eq. . For the 
classical-noise channel, this can be seen by expressing p' 
in terms of coherent states, i.e. 



Mx(|o)(o|) 



n + l \n+l 



(Al) 



(A2) 



The same relation applies for the thermal-noise channel, 
as can be seen by using the decomposition (1201) and the 
fact that the pure-loss channel maps the vacuum state 
into itself, viz., 

^(|0)(0|) =JV (1 _, )JV [£°(|0)<G|)] =^ (1 _ J?)iV (|0)(0|) . 

(A3) 



The entropy of p' is 



S(p' ) = -Tr 



Po a 



t^aln ■ 



M 



M + l 

(M + 1) ln(M + 1) — M In M , 



ln(M + 1) 



(A4) 



where we have used the fact that the average photon 
number of p' is M. 

APPENDIX B: PROOFS OF (|2SJ> AND (1271) 

In order to prove l|25|l and (|27|l . we employ the ampli- 
fying channel defined by the field transformation 



i/k a + — 1 



(Bl) 



where n ^ 1 is the amplifier gain and c is the annihilation 
operator of the amplifier's spontaneous-emission mode. 
With c in the vacuum state, the symmetrically-ordered 
characteristic function of the amplifying-channel map A K 
is easily shown to be 



(B2) 



Linking Eqs. (|17fl and l|B2|) . we find the decomposition 
rules 



cN _ cN' a 



for rj ^ r] , 
Af n = £°_ n o Aij^ n ) for n sC 1 , 
N n = A l/n o £ ° for n = (1 - rj)/^ 

where in I|B3J| . 

(1 - rj)N + rj' - rj 



N' 



1-T]' 



^ N 



(B3) 
(B4) 
(B5) 



(B6) 



Bound 1)25(1 now follows from combining Eq. I|B3(I with 
(|11[) . which establishes that minimum entropy is in- 
creased by concatenation of two maps. Using rela- 
tion (|B4|I together with Eqs. ijTSI) and lf2*U|) . we obtain 
the identity 



Nn = Mn,-n' ° £%- n > ° A/(l-n' 
P (n-n')/n' . 



(B7) 



which applies for n' £ [0, min(l, n)}. Bound (|27|l follows 
by removing the amplifier map using inequality i|ll|) . 



APPENDIX C: DERIVATION OF LOWER 
BOUND (l39t 



In this appendix we derive lower bound (|39|l for the 
minimum output entropy §(A< r ri ). This bound arises from 
the connection between the von Neumann entropy and 
the Renyi entropy of order two. Consider the family of 
states p with Tr(p 2 ) = t. As discussed in [2]J, the min- 
imum values of S{p) on this family are obtained from 
states that have a non-degenerate eigenvalue Ao and a 
fc-fold degenerate eigenvalue Ai = (1 — A )/fc ^ Ao, i.e., 



S(p) > -A o lnAo-(l-A )ln 
(1-Ao) 2 



1- A 



t 



A3 1 



k 



(CI) 
(C2) 
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Equation (|C2|) can be solved under the constraint 
Ai Aq with the result being 



A = 



1 - v /l-(fc + l)(l-H) 



for t e [l/(k + 1), 1/fc]. Substituting Eq. (C3j into the 
right-hand side of Eq. IjClfl . we obtain S(p) > F(t), with 
.F(t) being the function we have plotted in Fig. Ap- 
plying this result to the channel's output entropy, we find 
that 

S(N n (pj) > F(Tr[(AA„(p)) 2 ]) > F(l/(2n + l)) . (C4) 

The last inequality is derived by observing that F(t) is 
a decreasing function of t and that the minimum output 
Renyi entropy is achieved by a vacuum-state input , so 
that 



2 

N 



3 



4 

N 



k 



1 \ 








(C3) 




r"> ^ 





> 1 



FIG. 11: Beam splitter array needed to prove il-'il » . The in- 
put and output ports are numbered so that the jth input is 
facing the jth output: the input port 1 is fed with state p, 
the other k — 1 ports are fed with identical thermal states 
r of average photon number N. The transmissivities {r)j} 
are chosen in such a way that at each of the k output ports 
(indicated by arrows) we find the thermal-noise CP-map £^ 
with r\ — 1/fc. This corresponds to choosing rji = (k — l)/k, 

!J2 = (*- 2)/(* - 1), ..• ,»,»_! = 1/2. 



Tr[JV n (p) 2 } < l/(2n + l) 



(C5) 



Finally, follows from (JdJ, because F(l/(2n + 1)) 
coincides with the function on the right-hand side of (|39|1 . 



CD 




FIG. 10: Function F(t) defined in Eq. fCfy. 



APPENDIX D: DERIVATION OF LOWER 
BOUNDS (O AND (fTTTTi 

Here we derive lower bounds l|5U|) and l|51|) on the min- 
imum output entropy §(£^). 



Proof of lower bound (1501) 

We first prove that that (|50|l applies for 77 — 1/fc for 
all integers fc, and then we generalize to all 77. 

Consider the beam splitter array shown in Fig. 1111 
in which fc — 1 beam splitters of transmissivities 



rji , r/2 , • • • , rjk-i are connected in series and fed with fc — 1 
identical thermal states, each with average photon num- 
ber N. The transmissivity between the input port 1 and 
the jth output port is given by ej = (1 — r)j)r)j-i ■■ - rji. 
The beam splitters are chosen so that ej = 1/fc for all j, 
i.e., the jth beam splitter has r\j = (k — j)/(fc — j + 1). 
For example, with fc = 3 we have two beam splitters with 
transmissivities rji = 2/3 and 772 = 1/2, respectively, so 
that the transmissivity from the input port 1 to each out- 
put port is 6j = 1/3. The composition rule l|19|) can now 
be used to show that when the array is fed with a state 
p, at each of the output ports (apart from an irrelevant 
phase factor |2^|) we find the state £^ k (p). The output 
entropy of the joint state of all the outputs is equal to 
the total entropy of the fc inputs, because they are con- 
nected by a unitary transformation. For p a pure state, 
this entropy is given by the sum of the entropies of the 
thermal baths, i.e., (k-l)g(N). The subadditivity of the 
von Neumann entropy [Fj] implies that this quantity is 
less than the sum of the entropies of the single outputs, 



kS(£? /k (p))>(k-l)g(N) 



(Dl) 



which proves l|5l)|l for 77 = 1/fc. The case 77 ^ 1/fc fol- 
lows immediately by using rj' = 1/fc in (|24() and applying 
inequality HUll) . Lower bound (|5U|) for 77 1/fc is es- 
tablished by using rf = 1/fc in l|25|) and again applying 
inequality I|D1|) . 



Proof of lower bound (1511) 

As in the previous case, we first prove the bound for 
transmissivity 77 = 1/fc with integer fc and then we extend 
the proof to arbitrary 77. Consider now the beam split- 
ter array depicted in Fig. it is the scheme considered 
previously with classical noise maps added at each of the 
beam splitter outputs. The composition rules Eqs. H19|) 
and l|2U|) show that the output state of each port is the 



same |28|. and is given by £ ± 



N+nk/(k- 
/k 



1) 



(p). To calculate 
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FIG. 12: Beam splitter array needed to prove (1511 . Here 
the signals exiting each beam splitter encounter classical-noise 
channels. The transmissivities } are chosen as in Fig. 1111 



the joint state of all the outputs, we define a,j to be the 
annihilation operator of the jth input and V the unitary 
operator associated with the beam splitter array. This 
operator is given by U\U<x ■ ■ ■ Uk-i, where Uj is the jth 
beam-splitter operator, defined in Eq. JQ), which couples 
the mode Oj+i with one of the output modes of beam 
splitter j — 1. V produces the following field transforma- 
tion 



aV = a-A 



(D2) 



where A is the kx k real unitary matrix with |A,-j| being 
the effective transmissivity between the ith input and the 
jth output. The sign of Ay depends on the reflections en- 
countered by the field. The sign convention that follows 
from Eq. JQ) is that fields propagating from left to right 
(see Fig. 112ft acquire a 7r-rad phase shift when reflected, 
whereas fields propagating from top to bottom do not 
suffer any phase shift when reflected. This convention 
implies that Ay is negative if j > i and positive other- 
wise (e.g., A 23 = — v/(l - ?7i)(l - V2), A32 = 0). Given 
the input state R — p ® t ® ■ ■ ■ ® t (t being a thermal 
state with average photon number N), the joint output 
state of the Fig. 1121 map is obtained by acting on R first 
with V and then with the classical noise maps, i.e. 



R'= / d 2 pP(p) D(p)V RV^ D^p) 



(D3) 



where D(u) = exp(/2 • a ' — a ■ p^) with p = (pi, ■ ■ ■ , pk) 
a = (ai, • • • , dfc), and 



P( A /)=exp(-| / 2| 2 /n)/( 7 rn) 



(D4) 



Using Eq. l|D2p and performing a change of integration 
variables v = u ■ h.\ the output state can be written as 



R' = V 



J & 2 vP(v)D(v)RDHv) 
= V[N n ( P )®N n (T)®---®N n (T)]vK (D5) 
The entropy of this state is simply given by 

S{R') = S{M n {p)) + {k-l)g{n + N) , (D6) 



where we have used the fact that Af n (r) is a thermal 
state with n + N photons on average (see Sec. Ill h>|) . The 
subadditivity of the von Neumann entropy implies that 
S(R') cannot be greater than the sum of the entropies of 
the individual output states, i.e., 

k SiS^^ip)) > S{N n (p)) + {k-l) g(n + N) , 

(D7) 

which applies for any input state p. Note that when 
n = we recover inequality I|D1|) . as expected. Lower 
bound (|5 1 p - instead, derives by choosing N = 0, so that 



1 



g(n) . (D8) 



As in the previous subsection, we can apply the compo- 
sition rules (|24|l and l|25(l (using 7/ = 1/k) to extend the 
bound (|D8|I to any value of 77, obtaining inequality l|51[) . 



APPENDIX E: DERIVATION OF THE 
OUTPUT-STATE REPRESENTATION (TTfl) 

In this final appendix, we derive the Fock-state rep- 
resentation of the classical-noise channel's output state 
when its input is an arbitrary pure state. Using the Fock- 
state expansion of the input state, 



E 



4>n\n) 



(El) 



in Eq. I|17[l we immediately obtain the symmetric charac- 
teristic function at the output of the classical-noise chan- 
nel: 



X'(M) 



-n| M | 2 



jia' -ji a 



= e -(»+i/2)M' 

OO 




(E2) 
(E3) 



7^r fe (w 2 ) 



(E4) 



where {L^(z)} are the Laguerre polynomials. Recovering 
the output state p' from this characteristic function via 
J d 2 p x'{^) -C( — m)/ 71 "! an d performing the integration in 
polar coordinates, we obtain the desired Fock-state rep- 
resentation of p'\ 



(k + l\p'\k) = ^j 
(j + k + l)\ 



M ^ 



Cj + 



i\k\ (1 + n y+ k + l + 1 

x F(-j,-k;-(j + k + l);l- 



(E5) 



for k,l > 0, where F(a, /3;7; z) is the hypergeometric 
function. Note that (|E5|I becomes diagonal when the in- 
put is a Fock state. In this case, (|E5|I can be reduced 
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to H70|) by means of a transformation formula for the This work was funded by the ARDA, NRO, NSF, and by 
hypergeometric function and the connection between hy- ARO under a MURI program, 
pergeometric functions and the Jacobi polynomials |25j . 
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